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Abstract 

The non-perturbative behaviour of macroscopic loop amplitudes in the exactly 
solvable string theories based on the KdV hierarchies is considered. Loop equa- 
tions are presented for the real non-perturbative solutions living on the spectral 
half-line, allowed by the most general string equation [P, Q] — Q, where P 
generates scale transformations. In general the end of the half-line (the 'wall') 
is a non-perturbative parameter whose role is that of boundary cosmological 
constant. The properties are compared with the perturbative behaviour and 
solutions of [P,Q] = 1. Detailed arguments are given for the (2,2m — 1) mod- 
els while generalisation to the other (p, q) minimal models and c = 1 is briefly 
addressed. 
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1 Introduction 



Soon after the seminal works on non-perturbative 2D gravity appeared [EJ, a 
number of authors discussed the behaviour of macroscopic string amplitudes in 
those theories |^|, |[ || . In genus perturbation theory these are the surfaces with 
boundaries whose lengths scale in the continuum limit. It is possible however to 
discuss such correlators independently of this weak coupling expansion. Since 
the latter is only asymptotic and typically not Borel resummable one must input 
some non-perturbative information in order to even begin discussing such ques- 
tions. The mathematical possibilities are of course in principle endless. It is not 
so easy however to come up with a systematic formalism. In what follows the 
extra information that is required will be taken to be the principle that the KdV 
structure organising genus perturbation theory is present non-perturbativeljj^. 
This point of view was taken to its logical conclusion in refs.[f| |(| where it 
was found that new non-perturbative solutions were allowed that had not been 
previously considered. It is expedient to briefly review the main points of that 
work. Consider the (2, 2m— 1) minimal models, for which one has the KdV flows 
dt k u = — dt Rk+i[u] (notation is hopefully standard), which define the isospec- 
tral symmetries u — » u + edt k u leaving invariant the spectrum of the hamiltonian 
operator — 9 t 2 Q + u. Once the string susceptibility u has been determined these 
flows yield correlators of the local operators with couplings t^. The minimal 
physical requirement that one can impose on u is its scaling equation under a 
change of the length scale in the theory. For the set of parameters {u, {tk}} the 
KdV flows determine this equation to be ^ ; 

LK = (1) 

where 

11 = ^2(k + l/2)t k R k 

k>0 

= t ~j2( k + l / 2 w lL ^ lu ( 2 ) 

k>l 

D = V ± , L = -\D* + \uD + l -Du (3) 

and v is the renormalised value of 1/N. The string equation ([l]) may be written 
as the canonical commutation relation on the half-line, [P, Q] = Q, for co- 

1 This choice clearly has aesthetic and computational qualities but is nevertheless not, at 
present, a scientific argument. 
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ordinate (Lax) operator Q = D 2 — u and P = J2k>o tkQ^ 1 ^ 2 the generator of 
scale transformations of these co-ordinates, in a double-scaled Dyson gas with 
polynomial potential. 

To reproduce the ^-perturbative results as specified by the matrix models 
one must choose a solution which is asymptotically a solution of TZ — at 
to — ► — oo. One possibility is to choose an exact solution of TZ = 0, for which 
there are two types of real solution u. Those accumulating poles as to — + oo JjJ 
(type 1), and solutions of k odd critical points and non-singular flows around 
those points Q (type 2), for which u has a real asymptotic expansion at to — > oo. 
There is one other type of real solution to ([!]) with correct asymptotic behaviour 
as t — > — oo (type 3). This has no polesQ and has u — > as t Q — > oo (TZ — > to). 
It is realised by a Dyson gas restricted to lie on the spectral half-line i.e. there 
is an infinite potential 'wall' at the critical edge of the charge distribution. In 
the following section this representation will be used to derive Dyson-Schwinger 
loop equations. Figure 1 shows the numerical result for the type 3 susceptibility 
for pure gravity found in ref . Q . 

2 Type 3 Loop equations 

The scaling symmetry that led to (|l|) is one of the non-isospectral symmetries 
of the KdV hierarchy (see e.g.||). Together with the isospectral symmetries 
(KdV flows) it implies a family of constraint equations by applying the recursion 
operator LD^ 1 ; 

(LD- 1 ) n LTZ = , n>0 (4) 

These take the form of Virasoro constraints L n r = 0, u = —2D 2 log r ||To|| . 
There is an additional non-isospectral symmetry which, if applied, constrains the 
solution u further. Invariance under this transformation is the L_i constraint 
D1Z = 0, which is the hermitian matrix model string equation. This is the 
Galilean transformation: A translation of the Dyson gas potential, given in 
infinitesimal form by; 

u > u = u — e 

t k — ► t k =t k + e(k + 3/2)t k+ i (5) 
is an invariance of TZ but is not respected by ([!]) . 

2 The pole-free nature of the solutions was verified for the k = 1, 2 critical points in ref. fl| 
and has also been checked for k = 3 and the Ising model Q 
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In the case of type 3 there is a further parameter in the theory which has 
so far been neglected ||. For the solution on the half-line one should allow 
the infinite potential wall to be at some arbitrary scaled position, a say. This 
modifies the canonical position and momentum to Q + a and P + aP, where 
P = X)fc>i ^fc*2+ is the translation operator. In this way the string equation 
is modified to 0; 

[P + aP, Q + g] = Q + a (6) 

equivalently 

LIZ - <tDK = (7) 

Since a has the same dimension as u and (Q) is the scaling equation one identifies 
the evolution with respect to a: vd a u = —D1Z. This is the differentiated form of 
the L_i condition in the present case: L_it = d a r. The Virasoro constraints 
arc the expression of diffeomorphism invariance of the spectral line and the 
presence of a 'wall' has induced a boundary term on the right hand side. Similar 
boundary terms appear in the other constraints when a ^ 0, since varying the 
boundary as a — > a + ea n+1 implies that 

L„x = ,8) 

The Galilean and higher KdV symmetries are now respected when the correct 
transformation properties of a are taken into account. In particular to (|5|) one 
must add a — ► a — e. (^) becomes an equation for u in terms of t, independent 
of a. The constraints (^) may be rewritten in the old form for n > by using 
the invariance of the Virasoro algebra under 

L n — > L n = e-^-'Lne^- 1 (9) 

Given L_ir = d a r the equations (@) are easily seen to be equivalent to L n r = 0. 
The transformation (|]) describes a finite Galilean transformation by a, and the 
possibility of making redefinitions of this sort can be viewed as another reason 
for introducing a in the most general framework. 

The parameter a transforms in the same way as t m -\/t m under (|5|), where 
tfc = 0, k > m and t m is invariant in this case. Following the reasoning of ref. p| 
it suggests that its role is that of boundary cosmological constant. Indeed if one 
takes the usual expression for the renormalised macroscopic loop wavefunction 

|; 

poo 

<w(l)>= <z\cxpl(v 2 dz - u)\z> dz (10) 
J -t 
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a finite Galilean transformation by a exhibits the e~ al dependence of the loop. 
The net effect of the transformation is to define away the non-exponential de- 
pendence of < w(J) > upon a. In this case one can now say that a couples to 
the boundary operator jll| in the sense that 

^)_<n^)>=E^ < n^) > ( n ) 



where 



I), - (£).-£(»♦ vfcab m 



k>l 



da™ 



'nl 



(13) 



It is worth noting that there is a certain amount of ambiguity in the definition 
of the loop wavefunction, the choice ( |Io|) being perhaps the simplest. Two 
reasonable requirements might be that it agree with the results of Liouville 
theory and have a local (KdV) operator expansion at small /. The former is 
meaningful perturbatively in v and the latter perturbatively in I. This does not 
fix contributions non-perturbative in both of these parameters however. In this 
letter the definition ([l0|) will always be assumed. 

The arguments up to now have been somewhat heuristic so it is instructive 
to give a more careful treatment of some of the points. First a proof of the 
fact that (0) is the string equation of a Dyson gas on [a, oo) will be given. 
To agree with the conventions for non-universal constants adopted implicitly 
earlier, the following calculation is performed with an even polynomial action 
NV/A on the interval [—2,2], the ends of the charge distribution coinciding 
with infinite walls at ±2 (in un-scaled variables). The scaling regions around 
±2 thus furnish identical copies of the system on the half-line. Introducing an 
infinitesimal cutoff 5, let the scaled positions of the walls be ±2 =p cr5 2 . Further 
renormalised parameters are defined by; 

A = 1 + t Q 5 2m — = v5 2m+l 
An 

— = l-zS 2m Rn = l-u{z)S 2 , n~iV (14) 

for the neighbourhood of the mth critical point. The equations of motion in the 
orthogonal polynomial formalism |l2] ] are; 

(2n + l)A , 2A, 9 , s 

Jf <n \ XV \ n > = ^-^ 2 )P 2 (15) 
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— - v^<n-l|V \n> = — V^P n P„_i (16) 

The notation P„ is shorthand for P ra (2)exp(— NV(2)/2A) where P n (X) are poly- 
nomials orthogonal on [—2,2]. Terms involving P„ in (]l5|)(|l6|) are boundary 
terms picked up when one integrates by parts standard identities to derive the 
equations of motion. The left hand side of (|l^) is; 

(2n + l)A 
N 

2A 



- y/R n+1 <n\v'\n + l> -y/R^<n- 1\V \n> (17) 



y WRn+lPnPn+l + VKPnPn-l) (18) 

using (fL6|). From the work on polynomial potentials 0; 

nA ..,„/, ,2m. 



v <n-l\V \n>^R n = 8 2m (K(z)+0(5 2 )) (19) 
From ©@ and (£§) one has; 

6 2m Y 2 = ^(2-<j6 2 )P 2 (20) 

y = {n{z) + n(z + v 5) + o{8 2 )f/ 2 (21) 

and eliminating P„,P„_i from (|l8| ) gives the string equation at the first non- 
trivial order in 5; 

2 2 

^IZli' - V —{li) 2 -{u-o)K 2 =() (22) 

Differentiating once yields (0). The Dyson gas also supplies the appropriate 
boundary conditions u — > a as to — ► oo since u marks the edge of the charge 
distribution in the leading WKB approximation. In fact a BIPZ analysis [fli] , 
^| shows that in this approximation the charge density aquires a square-root 
divergence at the wall (t > 0). More generally at to < this divergence has 
exponentially small residue, as can be most easily seen from the conventional 
form of the Dyson-Schwinger equations. These can be derived in the usual way 
[ fl3"| since the type 3 solutions have a path integral representation, and should 
correspond to the Virasoro constraints described earlier when Taylor expanded 
along the lines of refs.[flO|. Starting from the partition function on [— S, E] say, 
in unsealed variables; 

Z= HdX, A 2 (A)e- JV S. y(A ' ) (23) 

J — S - i 
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and denning the loop generating function and its Laplace transform (marked 
loops are used for convenience); 

i 

#<n = ^rh-, < 25 > 

i 

one may perform an infinitesimal change of variables A; — * Aj+e/(T— A^) in (p3|). 
The only new contribution to the standard analysis |T^| is from the variation of 
the boundaries, giving 

dZ ( I 1 \ 

e— I - =— = I (26) 



<9£ \T -E T 
The first loop equation is then 

/ 2£ dZ 

V (T) <W(T)> +H(T) =< W(T) W(T) > + 2 _ m (27) 

where H(T) is a linear combination of < W(0) >,< W (0) >, .... In terms of 
W(L) one has the equation of motion (introducing connected correlators); 

V '(jfc) <W(L)> C = £ dL' (-L<W(L')W(L-L')> C 

+ <W(L') > c < W(L - if) >„) - 2si ^ S fl %fe ) 

showing a new source-like term for the loop generating function. 

Following ref.jj] it is a simple matter to take the continuum limit by intro- 
ducing, in the neighbourhood of the mth critical point, renormalised parameters; 

S = S c — a5 2 , T — -S c - t5 2 

<W(T)> C = hy'i^ + S 2 ™- 1 <w(t)> c (29) 

M M 
<n^m)>c = ^ m+ 2-M( 2m+ 3) <JT^ (Tj)>c 

1=1 1=1 
Equation ( p7|) becomes 

, , r> 9 , x , n / \ v 2 d log Z , . 

< w(t) > 2 C +v 2 < w{t)w(t) > c = tt(t) + ^— (30) 

a + r a<7 

For completeness the higher loop equations are 

M M 

2<w(t) > c <w(t) J[ w(n) > c +v 2 < w(t)w(t) J[ w(n) > c 

i=l i=l 
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+ <w(T)Y[w(T l )> c <w(T)Y[w(T :j )> c 

I ,J iel jej 

E~ / x 9 w(Ti) — w(t) _ . . 
<» ri . . . ^ ^ . . . w(t m ) > 
riT- T: — T 



<O a \\w{n)> c +^—^- < Y[w{t 1 )> c (31) 



M 9 g M 

, . da 

i=l i=l 



where Oq is the (bulk) puncture operator and the index sets /, J are such that 
I U J = {1, . . . , M}: I, J 7^ = I n J. 7t(t) is a polynomial in r, of degree 
2m — 1 at the mth critical point, which determines the structure of the charge 
density in the scaling region, given by ^Disc < w[r) >„=o in the saddle-point 
approximation. The last term in ( |30| ) should be ^r^r / TZdto from the discussion 
following (0), which is exponentially small in v at to — oo but is 0(1) for 
to — > oo. The physical meaning of the divergence as r + a — > + , is that 
large loops become unsuppressed, since w(r) is the Laplace transform w.r.t. r 
of the renormalised macroscopic loop wavefunction (|To|) . Because of the e~ al 
dependence of (|l^) the Laplace transform only exists for r > — As indicated 
in the discussion following ( |2^ ) , a concomitant divergence appears in the charge 
density as r + er — > 0~ . 



3 Discussion 



The previous calculations have shown that type 3 solutions yield very simple 
e~~ al behaviour for macroscopic loops. This exponential behaviour is much like 
the ^-perturbative result. At each order of a WKB expansion the charge density 
has support on the half-line (— oo, —u] of the real r axis, specified in the leading 
approximation by the (single) cut in ^tt(t). The rest of the real r axis describes 
the loop function w(t) which on inverse Laplace transformation is seen to have 
e~ ul times power law behaviour. At genus zero w(l) contains universal terms 
with inverse powers of I, as emphasised in refs. |fi"5| , since 7r(r) is polynomial: 

< w (l) >= - klt k r k -^ 2 + 0(l 1/2 ) (32) 

v fc>0 

The Laplace transform does not exist because of these terms, but one can pro- 
ceed by differentiating w.r.t. r a sufficient number of times. It is in this sense 

3 In fact there are also negative powers of I in w(l) preventing naive Laplace transformation, 
but these come from low genus and can be systematically isolated. 
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that w(l) and w(t) are transforms of one another. Alternatively one can work 
at fixed 'area' instead of fixed bulk cosmological constant. The offending terms 
are analytic in the latter and so do not contribute to an inverse transform to 
fixed 'area'. They correspond to finite loops spanned by infinitesimal surfaces. 

The z/-non-perturbative exponential behaviour in I of loops for type 1 and 
type 2 solutions is more complicated. For type f the discreteness of the spectrum 
of — D 2 + u implies that < w(l) > behaves like an infinite sum of exponentials 
with different arguments Q . Only in the / — » oo limit does one recover a simple 
e -eo( behaviour, where eo is the lowest eigenvalue. It is sometimes suggested 
that a solution satisfying a loop equation such as (^0|), derived from a path 
integral representation, is 'physical' and that, by implication, one that does not 
is 'unphysical'. There is no known path integral formula for type 1 solutions. 
This is not (presently) known to contradict any physical principle however. The 
hcrmitian matrix model is the path integral representation of type 2 solutions, 
and shows that the loop expectation is always diverging as / — > oo for such 
solutions H . This markedly different behaviour is due to the fact that the charge 
density has support on the whole spectral line, in particular it has an exponential 
tail. For a tail p(e) ~ exp(— |e| p ) as e — » — oo a dimensional argument shows 
that <w(l) >~ expZ 1+1 /( p_1 ) as I — > oo. This means that one cannot define the 
loop function w(t) in this case. 

To conclude one may note some possible generalisations of the results of 
this letter to other models with c < 1. For the (p, q) minimal models described 
by the generalised KdV hierarchy the string equation (scaling equation) 
[P, Q] = Q will provide new solutions analoging those of type 3. Generally one 
only knows how to treat macroscopic loops embedded at a single point in the line 
of q— 1 points. By the same argument [ pd| , shifting the non-derivative part of Q, 
one can identify a parameter coupling to a boundary operator. Perturbatively it 
appears that only one parameter can be generated in this way (e.g. for the Ising 
model (4,3) it is the boundary magnetic field), which led the authors of ref. [pT| 
to conclude that certain operators could not be expressed in the KdV formalism. 
However in the case of type 3, non-isospectral evolution equations also play a role 
in defining couplings and one can imagine restricting the q — 1 types of charge 
to different half-lines i.e. the wall becomes 'time'-dependent. The necessary 
argument in terms of the W-constraints to confirm or deny the validity of this 
naive picture is a little involved. A detailed account of the [P, Q] = Q version 
of other minimal models will be given by the authors of ref. M . At c = 1 the 
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picture is similar. By placing a wall in the scaling region given by the inverse 
quadratic potential jlTj one has a stable quantum mechanical system providing 
a non-perturbative definition of the theory. More particularly the macroscopic 
loop amplitude at time t is non-perturbatively well-defined; 



where the wavefunctions \& vanish at the wall A = a and the matrix model 
corresponds to the continuation to euclidean time. It exhibits a dependence 
cxp — la, at the expense of introducing a linear term in the potential. The 
extra parameter a more generally has a continuous argument since again it may 
be time-dependent. Understanding of the possible flow structure at c = 1 is 
still hazy. Is stabilisation by a simple wall, possibly fluctuating in position, the 
only quantum mechanical problem non-perturbatively compatible with a flow 
structure organising perturbation theory, akin to c < 1? 
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Note Added: After this letter was typed a preprint appeared where 
quantum mechanics on the half-line is discussed with regard to the requirement 
of unitarity of tachyon scattering at c = 1. 
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Figure Caption 



Figure 1 : This shows the (unique) numerical solution of type 3 for pure gravity. 
The Hamiltonian —D 2 + u has continuous spectrum down to the value of the 
left asymptote of the potential u (zero in the figure shown). Those eigenvalues 
are more-or-less directly related to the positions of the Dyson gas charges on the 
spectral line, which should lie on the positive halfline. Note that although there 
is a small well in u, the previous identification indicates that it is too shallow to 
support bound states below the continuum. A rough numerical estimate using 
the proportions of the figure confirms this. As is explained in section 2, more 
generally the left asymptote can be u = a if the Dyson gas is restricted to 
[er, oo). It is possible that non-perturbative subleties arise as a turns negative 
e.g. through an instability in u similar to that found in ref. @). 
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